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Abstract 

The representation of a magnetic field by the cross-product of the gradients 
of two scalars has recently seen wide use in plasma physics and in the study of 
energetic particles in space. The properties of such a representation are re- 
viewed, and as an example of its application, the first-order guiding center mo- 
tion of a charged particle in a time -independent magnetic field is derived in 
canonical form. 
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"From my earliest experiments on the relation 
of electricity and magnetism, I have had to think and 
speak of lines of magnetic force as representations 
of the magnetic power; not merely in the points of 
quality and direction, but also in quantity" 

Michael Faraday, Experimental Researches 
in Electricity, § 34. 

"As I proceeded with the study of Faraday, I 
perceived that his method of conceiving the 
phenomena was also a mathematical one, though 
not exhibited in the conventional form of mathe- 
matical symbols. I also found that these methods 
were capable of being expressed in the ordinary 
mathematical forms..." 

James Clerk Maxwell, Preface to "A Treatise 
on Electricity and Magnetism", 1873 


I. General Properties 
Introduction 

The concept of magnetic field lines is to a large degree due to Michael 
Faraday. Faraday's views were often more intuitive than mathematical, and he 
regarded magnetic field lines (or to use his term, magnetic lines of force) as a 
useful way of visualizing the magnetic field. Later generations of scientists, 
more inclined towards mathematical abstraction, replaced Faraday's qualitative 
description with the more precise terms of vector fields and of vector and scalar 
potentials, relegating the use of field lines mainly to elementary textbooks. 

Recently, however, there has arisen renewed interest in such lines, following 
investigations of charged particle motion in magnetic fields and of plasma 
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phenomena, effects in which the configuration of field lines is significant. The 
purpose of this article is to describe a mathematical tool which is appropriate 
for such applications and to review some of its uses. 

Euler Potentials 

A general vector field in three dimensions requires 3 scalar functions of 
position for its description, but the magnetic field B requires only two, since it 
satisfies V • B = o (in this work, B rather than H will be regarded as the field 
vector). One such representation (which can be shown to be generally possible 
at least locally) is 


B = Va x V/5 (1) 

The scalars a and /3 are termed Euler Potentials. They naturally lead to 
the vector potential 


A = aV/3 

satisfying the gauge condition (A • B) = 0. Clearly, they are far from unique, 
for an arbitrary function of a can always be added to /3, or vice versa. More 
generally, given the Euler potentials (EP) a and /3, one can show by inspection 
that any pair of functions u (a, /3) and v(a, J3) may replace them in Equation (1), 
provided 


3 (u, v)/B (a, /3) - 1 (2) 

A necessary (though not sufficient) condition for a and /3 to satisfy eq. (1) 
is that they be independent solutions of the same linear partial differential equation 
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B • Vo. = 0 

B ’ Vyfl = 0 

Since Va and V/3 are perpendicular to B, surfaces of constant a and /3 are 
tangential to 'he field at all points, and the same hold j for the lines along which 
such surfaces intersect. Such lines are therefore magnetic field lines. 

Suppose that we know the EP representing a given field in some region in 
space. We then have two families of surfaces 

a(x, y, z) = a. 

/3(x, y, z) = /3. 

Each field line in the region is the intersection of two surfaces, one from 
each family. It is consequently characterized by two parameters (a. ,/3. ), 
equaling the constant values assumed by a and (2 along it. In this manner, a 
formulation of the field in terms of EP affords a direct representation of the 
field lines, in a way not possible with the ordinary vector potential A. 

A Short History 

Long before the nature of magnetic fields was understood, mathematicians 
investigated the velocity field v of incompressible fluids , which is likewise 
solenoidal, since the equation of continuity then reduces to 

V • v = 0 

The field lines of v are usually termed streamlines. Leonhard Euler (l) 
was the first to introduce into the description of v "stream functions" F and G 
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which are conserved along streamlines (Figure 1). Because of the difficulty in 
deriving such functions explicitly, subsequent developments of fluid dynamics 
showed relatively little interest in the representation of streamlines. An excep- 
tion was the Stokes stream function (2) , derived by Sir George G. Stokes early in 
his scientific career <3) , which gives streamlines in the more tractable case of 
axisymmetrical flow. 

The first application of EP to magnetic fields is due to Sweet (4) and was 
continued by Dungey (5) . In a work on solar magneto-hydrodynamics, Sweet 
represented the field as 

B - F(V<£ x Vi/,) (4) 

with F a function of position. The functions <p and are evidently conserved 
along field lines, and must therefore be functions of the EP. Indeed, if one re- 
places (a, /3) in eq. (1) by two functions 4>(a, /3) and ^(a,/3) not bound by eq. (2), 
one obtains a relation of the above form, with F equaling the jacobian. Pairs of 
functions such as <p and 4> will be termed unmatched EP: the cross product of 
their gradients is parallel to B, but not proportional to it in magnitude. 

The use of EP as given in eq. (1) originated in investigations of charged 
particle motion by Northrop and Teller (6) , Grad (7,8) , Gardner Ray (10) and 
others. It has become customary to designate them in that case by a and /3; if 
one of the potentials can be derived in a straightforward manner from considera- 
tions of symmetry (e.g. the azimuth angle <P in axisymmetrical fields) the 
symbol /3 is usually reserved for it. 
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While there has been considerable uniformity in notation, less of it exists 
in terminology (and on one occasion, at least (11) , the symbol a has been used as a 
name). Grad has referred to a and as "Euler Stream Functions" (7) and "Flux 
Coordinates" (8) , while Truesdell, in a work on fluid dynamics (12) , has called them 
"Euler's Potentials". Most authors, by and large, have studiously avoided using 
any special name. 

Mention may alsj be made here to the term "Monge Potentials", applied to 
functions (a ,0,y) in the decomposition of an arbitrary vector field V in the form 

V = (Va X V/3) + Vy (5) 

which is implied in a 1784 work by Gaspard Monge. By this definition, a and /3 
could be called Monge potentials. However, the application of eq. (5) to divergence- 
free fields is particularly ambiguous. For instance, the main geomagnetic field 
created by sources in the earth's interior can be represented either as 

B = - V y (6) 

(the customary representation, with y expanded in spherical harmonics) , or by 
means of a cross product as in eq. (1), and of course by a large variety of ex- 
pressions containing terms of both types as well. 

In order to avoid such ambiguities , it appears best to use a special term 
for a and (3 in divergence -free fields. The term Euler potentials, recently in- 
troduced in work on the geomagnetic field, seems to be the most appropriate 
here, giving credit to the originator of the formulation and at the same time 
stressing that it provides a representation of the field by means of auxiliary 
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functions, equivalent to that provided by the vector potential A or (in source-free 
regions) by the scalar potential y . 

Properties of Euler Potentials 

Given the magnetic field B in the region surrounding a point P , it is always 
possible, in principle, to derive a set of Euler potentials describing it in some 
vicinity of the point. The proof of this ( 13) will now be sketched. 

The linear homogeneous partial differential equation (3) admits, in the 
vicinity of P, two independent solutions. Let u and v denote two such solutions 
(obviously, any well-behaved function f(u, v) is also a solution). Since they are 
independent, they may be supplemented by a 3rd function w so that (u, v, w) can 
be used as curvilinear coordinates in the vicinity of P. 

Since V u and V v satisfy (3), their cross product is tangential to B and 
satisfies (4) 

B = F(Vu x V v ) 

The divergence-free character of B limits the choice of F: 

V • B - 0 r (3F/3w) Vw • (Vu * Vv) 

= (9F/ <3w) o (w, u, v)/B (x, y, z) 

Since the functions are independent, the derivative of F must vanish, so 
that F is a function of (u. v) alone. One may now form 

X'(u, v) = f F(u, v')dv‘ 

*^ v o 
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with v Q some arbitrarily chosen constant. Then 

V X - \Vu+F(u,v)Vv 

with A. some irrelevant function related to v 0 . The function x constructed in 
this manner is thus indeed the EP matching u, which completes the proof. 

For certain symmetrical fields, Euler potentials are readily obtained. For 
instance, for a poloidal axisymmetric field (e.g. that of a dipole), the vector po- 
tential may be chosen in the azimuthal direction 

A = A(r, tf) <f> 

leading to the choice 

a - A(r,6 l )r sin a 

- <P 

Similarly, one may represent by Euler potentials pure toroidal fields 

B - V x t\p - x V (r 2 / 2) 

and two-dimensional fields (with neither dependence on cartesian z nor a com- 
ponent in its direction) 

B - V x Az - VA(x, y) x Vz 

In more general cases, however, the analytic derivation of EP is quite diffi- 
cult, since they enter eq. (1) in a nonlinear fashion involving products of deriva- 
tives. This non-linearity prevents the superposition of solution (and thus the 
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knowledge of EP for Bj and B 2 does not help one find those of [Bj + B 2 ] ), except 
for the ''ase when the superposed fields all have one Euler potential in common 
— e.g., when they all belong to one of the symmetric classes previously listed. 

In fields deviating only slightly from symmetry, perturbation techniques may be 
employed, and these have been successfully used in connection with the repre- 
sentation of the main geomagnetic field (14) . 

The non-linearity is a distinct disadvantage compared with the ordinary 
vector potential A which may be derived by linear superposition of contributions 
from the various field sources, thus providing a general "brute force" method of 
numerical derivation. However, as will be shown in what follows (and may, in- 
deed, be inferred from the existence proof), somewhat similar methods also do 
exist for Euler potentials. 

Suppose that one is given the line pattern of a magnetic field in a certain 
region, as well as the field B on some surface a through which each field line 
in the region passes exactly once. In that case, we shall now show, there exists 
a numerical procedure for deriving a and which, indeed, could be the most 
straightforward way of obtaining B from such data. 

The problem, incidentally is not without practical interest. Consider the 
earth's magnetic field: one can observe it with great precision over the earth's 
surface, while satellite observations^ S) provide us with a general field line 
pattern (influenced by factors that are not too well understood) far away from 
earth (the satellites also give information about the field's magnitude, but due to 
the steep intensity fall-off and the appreciable variations observed, such informa- 
tion tends to be less useful). The problem of reconstructing the field on the basis 
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of such information is very similar to the one previously described, except that 
due allowance must be made for the fact that most field lines cross the earth's 
surface twice. 

Let (K, fjL, v) be curvilinear coordinates, such that the surface o is charac- 
terized by a constant value v 0 of v (for the case of the magnetic field observed 
outside a spherical earth, v will be the radial distance r or some function of it). 
Then surface points may be labeled by two functions of the remaining coordinates , 
e-g. by 

u = u(\, m) 

v = v(\,/x) 

Since the field line pattern is given, one may associate with each point in it 
two quantities U and V, equal to the values of u and v at the intersection with a 
of the field line passing through the point. By this definition, U and V are spatial 
functions, numerically derivable and having the property 

v o) ~ u(^,/x) 

„) = v(A fi) 

Both U and V are conserved along field lines and thus constitute a set of EP, in 
general unmatched. For some particular choices of (u, v), however, they will be 
matched, and it will be now shown that such choices can be obtained using only 
information about the field on cr. 

Given an unmatched set (U, V), let [U, /3(U, V)] be a matched one; then the 
grid on a- generated by [u, /3{u,v)] g< ne rates matched potentials. On a, every 


vector quantity may be resolved into components tangential and normal to 
the surface, denoted by subscripts m and i respectively, with the latter com- 
ponent oriented along V v. One gets 

B, = (B/S/BV) [VU M x VV„ x 

The gradient components involved in the cross product are readily derived 
from the known functions u and v, since (K , /x , v not necessarily orthogonal) 

vu„ - (BU/BP^ V\„ + (SU/3^ V Mn 

- (Bu/3A.) V\ m + (3u/3u) Vyu M 

and similarly for V. Since B on a is assumed given, one may derive a function 
of u and v (to which U , V reduce on a) 

f(u,v) = 

and obtain /3 by integration 

/3(u, v) = J* f (u, v 1 ) dv' 

This gives J3 within an arbitrary function of u , but such a function does not con- 
tribute to the cross product. 

The preceding existence proof and construction methods, unfortunately, are 
only valid locally, enabling one to derive (a, /3) only in some neighborhood of a 
given point and not necessarily for all space or even for the entire region of in- 
terest. When one attempts to extend such construction over larger regions, one 
often finds that the EP are no longer single valued. 
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As an example (16) , consider the field of a current in a ring filament. This 
field is axisymmetrical , its a -surfaces are toroidal rings nested inside each 
other (all of them enclosing the conductor) and its /3 -surfaces are meridional 
planes. By adding an infinitely long current filament perpendicular to the ring 
and passing through its center, the field lines acquire a slant, causing each to 
spiral around the ring (figure 2). 

Let a point P be given in this field, and let Euler potentials ( a , j3) be derived 
(numerically or otherwise) in some small region r surrounding it. This labels 
each field line with a pair of numbers, and one may extend this labeling to points 
outside r by simply following labeled field lines after they have left the region. 
Ultimately, however, these field lines will have circled the central wire and will 
have returned into r, encountering previously labeled points, in general with 
different values of (a, j3). Such points will then have raore than one pair of values 
(a, /S) associated with them; the number of such pairs can usually be made arbi- 
trarily large by taking into account field lines that have circled the wire more 
than once. 

Mathematically, this ambiguity can be sidestepped by introducing a cut — 
a surface of discontinuity at which the labeling is artificially terminated. In the 
example cited, the ^'ane surface (or any other) bounded by a suitable section of 
the ring and two radii can serve this purpose , so that each (a , fi) line begins on 
the surface and ends there as well. 

In actual practice, this artifice has limited value, for given problems often 
involve following field lines for long distances (17) (18) . For example, particles 
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trapped in a stellarator field (which resembles that of Fig. 2) may, under 
suitable conditions , execute many circuits around the central wire while following 
a field line. One may still use a surface of discontinuity to achieve unique char- 
acterization; on such a surface, then, to each (a, /3) line-segment corresponds 
some point P. After one circuit, the representative point is usually shifted some 
distance, so that one can regard each circuit as a mapping of one point on the 
surface into another, and repeated circuits as iterated mappings. The study of 
long-term behavior of field lines — at least, in situations where the above- 
mentioned shift is not too wild — is then transformed into an analysis of iterated 
mappings (19)(20) . 


As was noted before, Euler potentials are useful in problems where the 
physical situation involves the field line configuration. In such cases it often 
turns out that the equations are easier to handle when of the 3 coordinates of 
position needed to define a point in space, two are chosen as its Euler potentials 
a and fi. The 3rd coordinate is usually chosen as the distance s along a field 
line, measured from some arbitrary surface. 

By its definition, Vs has a component of magnitude unity along B. However, 
it is not generally parallel to B ; if it were , this would mean 

B = A. Vs 
B • (V x B) = 0 

a condition not always fulfilled. A transformation of (a, /3) will not in general 
affect s , since that latter quantity depends only on the field line structure, not 
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on its representation by Euler potentials. Similarly, without affecting the EP, s 
may be transformed by changing the surface from which it is measured. The 
general form of such a transformation is 

s - s ' = s + f (a, / 3 ) 

with f an arbitrary function. 

With these preliminaries, we now proceed to an example of the application 
of Euler potentials. 
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II. Application 


Guiding Center Motion 

As an example of the application of Euler potentials we shall now carry out 
the canonical derivation of the nonrelativistic guiding center motion of a charged 
particle, following a course outlined in general form by Gardner (9) . 


Consider a particle of change e and mass m in a time-independent magnetic 
field that is given in terms of Euler potentials , with canonical variables of posi- 
tion x. and their conjugate cononical momenta II. . The Hamiltonian will be 


H 


II - (e/c) 




/ 2m 


( 8 ) 


We assume that the Euler potentials a and /3, as well as the arc length s 
along field lines are known functions of the x. . Following Gardner, we introduce 
a canonical transformation to variables (Q A ,P i ) , generated by 


F(x.,P.) = sP 2 + /3Pj + aP 3 - (c/e) P, P 3 


( 9 ) 


From the transformation equations 


One obtains 


a 

0 


Q. = 3F/3P. 

= (c/e) P, + Q 3 

= Qj + (c/e) P 3 
s = Q 2 


( 10 ) 


(ID 

( 12 ) 

( 13 ) 
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while the remaining equations 


n. - SF/Bx, ( 14 ) 

yield, when substituted in (8) 

H = (l/2m) |p 2 V s - (e/c) Q 3 V/3 + P 3 V a | 2 (15) 

Homogeneous Field 

Assume at first that the field is in the z direction and has a constant in- 
tensity B. We may then take 

a = x 

P = By (16) 

s = z 

The Hamiltonian separates into two parts, representing motion parallel and 

perpendicular to the field 


H„ = (1/ 2m) Pj = E„ 

(17) 

*■/ 2m) | (eB/c) 2 Q 2 + P 3 2 } = E, 

(18) 


The parallel part simply gives constant motion along the z axis 

P 2 = m v t| = constant (19) 

while H, resembles the Hamiltonian of the harmonic oscillator and yields 

Q 3 = (2 E, /mco 2 )* sina>t = R g sin cot (20) 
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( 2 m E, cos cot 


(21) 


where 


P 


3 


co - e B/m c 


( 22 ) 


is the well-known gyration frequency and R g is the particle's gyration radius. 
The motion is periodic, and the associated action variable is 


J - fr, dQ 3 


277E, /co 


Let 


(23) 


W - W 0 (Q 3 ,J) - E t 


(24) 


be the solution of the Hamilton-Jacobi equation for H; W Q generates a transfor- 
mation to action-angle variables (J , ft) and satisfies 

(eB/c) 2 Q 3 2 + (BW 0 /BQ 3 ) 2 = 2mE, (25) 

from which (J substituted for E^) 

W Q = f |(JBe/7Tc) - (eB/c) 2 Q 3 2 }* dQ 3 (26) 


Differentiation beneath the integral sign gives 

ft = dvi 0 /dj = (277)' 1 arcsin Q 3 (rreB/Jc) H (27) 

which, compared to (20), identifies 277 ft as the accumulated gyration angle, so 
that ft grows linearly in time and increases by 1 each period, as required of a 
canonical angle variable. 
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Both Pj and Q t are absent from either of the Hamiltonians, showing that 
they are constant. By equations (11) and (12), they may be identified with (e/c) a 
and /? of the field line around which the particle spirals. In what follows, the 
point 

(a. As) r - ^(c/e)P,, Q r Q 2 j 

around which the particle instantaneously revolves will be termed its guiding 
center and quantities associated with it will be distinguished by subscript zero. 

Inhomogeneous Field 

In an inhomogeneous field, one can define u> and R g at every point by using 
the local value of B. Suppose next that the ratio of R g to the scale L on which the 
field varies — a ratio which will be loosely designated by e — is everywhere 
much less than unity; this occurs when either the extent of the inhomogeneity or 
the perpendicular energy E ( is small. Then to lowest degree of approximation, 
the particle on its local scale senses a homogeneous field, and its motion is the 
same as derived in the preceding section, with its canonical variables and 
guiding center defined as before. 

To obtain a better approximation, one has to take into account the variation 
of field quantities over the course of a gyration, and this is usually done by ex- 
panding them around the guiding center. If /? is the value of /3 at the guiding 
center, the expansion of a function f(/3) would be, by eq. 11 

f 0 ) = f(/3 g ) + (c/e) P 3 (3f/3/S) g + ... (28) 


17 


Here P 3 is of order R g while f /(Bf/B/3) is of order L, so that the ratio of 
the terms is 0(e) and, under the assumption made, the series provides a good 
approximation. The approximation is even better (by one order of e ) when we 
consider quantities averaged over one gyration, since in homogeneous fields 
such averages vanish and in slightly inhomogeneous ones it seems logical (and 
may be shown) that they are no larger than 0(e ). The gyration-averaged ratio 
of the terms in (28) is thus 0(e 2 ). 

By the same reasoning, since in a homogeneous field a and f3 of the guiding 
center are conserved , in a slightly inhomogeneous field their deviation from 
their initial values (a Q , /3 0 ) will be 0(e ), as long as the time scale is not too 
large. All this will be assumed a-priori in the following calculation but will 
also be borne out by the final results. 

To extend the canonical formalism to slightly inhomogeneous fields, one 
replaces the transformation generator of eq. (9) by 

F ( x i - p i ) = s o p 2 + /3p i + aP 3 " H ^ p 3 2 ” (c/e) P x P 3 (29) 

where s is replaced by 

s 0 = s - (a - a 0 ) a - (/3 - /?„) b (30) 

and 

a - B- 2 j(V/3) 2 (V s • V tt ) - (Va • V/3)(V S • V/3) j 

b B- 2 | (Va) 2 (V s • V/3) - (Va • V/3)(V S ■ Va)J (31) 

A. - (c/ e) (Va • V/3)/(V/3) 2 
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Retracing the transformation equations 


a ^ (c/e)P 1 <■ Q 3 + aP 3 (32) 

* 

P ~ Q, + (c/e) P 3 (33) 

s o " <? 2 (34) 

and 

H ( 1 / 2m) | P 2 Vs 0 - (e/c) (Q 3 * \P 3 ) Vp^ + P 3 Va - % P 3 2 VA J 2 (35) 

Le* H 0 denote that part of the Hamiltonian that is contributed by the scalar 
products (Vs 0 ) 2 , (Va) 2 , (V/3) 2 and (Va • V/5 ). Using 

B 2 = (7a) 2 (V^y 2 - (Va • V/l) 2 (36) 

one finds 

H 0 - (1/ 2m) j P 2 2 (Vs 0 ) 2 + (e/c) 2 (V/?) 2 Q 3 2 + P 3 2 B 2 /(7/i) 2 j (37) 

The remaining terms in H can be classified by the scalar products which 
they contain, of which there exist 6. It will now be shown that these terms are 
all of a lesser order of magnitude than those contained in H 0 . 

Consider first the terms involving (Vs 0 • Va) and (Vs 0 • V/3). if b is the 
unit vector in the field's direction, one finds 
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( 38 ) 


V Sq - B - (a - a 0 ) V a - (/2 - /3 0 ) Vb 
= B + 0(e) 

Thus these terms are 0(e) — and, since they turn out to have factors of Q 3 
or P 3 , their averages are 0(e 2 ). The reason for replacing s with s 0 now 
becomes evident: in a slightly inhomogeneous field, the term P 2 Vs in (15) 
usually does not represent the momentum component parallel to B , whereas 
the term P 2 V s 0 in (35) does, to the lowest order of approximation. 

The other terms all involve which is of order L" 1 , and they therefore 
belong to the same order, except for the term involving (V\) 2 , which is 0(e 2 ) 
even without any averaging. To this collection of 0(e) terms with 0(e 2 ) averages 
we new add several others by expanding the functions of position contained ir H 0 
around the guiding center in the general manner of eq. (28), using the relatr -.s 
(32) and (33). Only the lowest order terms of this expansion are retained in H 0 . 

If we wish to solve the motion in Q 3 and P 3 to order e — for instance, to 
derive the first-order correction to the magnetic moment (see later) — the 
higher order terms have to be retained and no averaging is permitted. Such a 
calculation, obviously, is quite lengthy. Here, instead, we shall only investigate 
the average first-order motion of the guiding center coordinates, in which case 
H 0 is all that is needed. We begin with a transformation to new canonical 
variables 

(J > pj> q r Pj> q 2 ) 

with a generator that is a cross between the one of eq. (26) and that of the 
identity transformation, namely 
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W(p.Q) = e(V/3)V cB j j^JBc/77e(V/3) 2 j - Q| J* dQ 3 + Pl Q, + p 2 Q 2 (39) 

where B and (V/3) 2 are defined for mixed guiding center variables 

a = cpj/e 

0 ~~ Qi 

s o = Q 2 

Applying 

P 3 = 3W/3Q 3 

one finds that the terms in H 0 involving P 3 and 0 3 condense to 

JBe/77c = ,uB (40) 

The expression obtained for fi resembles that of (27), except that B there 
is replaced by (V/3) 2 /B. 

To obtain the rest of the variables , one first has to derive the integral 1 1 
in (39). One gets two terms — one containing a factor Q 3 , the other proportional 
to ft. Fortunately, in the coefficient multiplying the latter term, the dependence 
on position cancels out, so that when one derives quantities such as 

P, = BW/BQj 

1 3 contributes no terms proportional to ft, only such ones that resemble in form 
the second part of (28), i.e. of order e with 0(e 2 ) averages (a possible vanishing 
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denominator also cancels). Apart from such terms, which will be neglected, 
the remaining new canonical variables are the same as the corresponding old 
ones. We get, to our order of approximation 

H - (l/2m) |p 2 2 (Vs 0 ) 2 + M B j (41) 

Since ft appears only in 0(e) terms, / j . is constant to the lowest order: this is 
the well-known adiabatic invariance of the magnetic moment to its lowest order. 
If p, is the perpendicular component of the momentum, then 

M = Pj*/B 

We now use (41) to derive the first-order variation of q, and pj, the so- 
called guiding center drifts. The velocity component normal to B is defined 

VjL = B~ 2 x ( v x B)j 

Using the identification of (cp/e, q x ) as (a ,/3) of the guiding center, Hamilton T s 
equations give 

(v x B) - Va(v * V/3) “ V/3( v * Va) 

- /3 Va - a V/3 

(c/e) (v x B) = Va(3H/Ba) + V/3(3H/B/3) 

= VH - Vs 0 (BH/3s 0 ) (42) 


with (38), this gives to lowest order 

v, te (c/eB 2 ) (B x VH) (43) 
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Referring to (41), we find that v, consists of two parts. The second term in H 
contributes the so-called gradient drift, proportional to VB 

v g = ( / ac/2meB 2 ) (B x VB) (44) 

The first term contributes the curvature drift, proportional to the curvature of 
field lines 


V c - (c/meB 2 ) p 2 |b x (Vs 0 • V Vs Q )j 
From (38), to lowest orders 

Vs Q • VVs 0 = B • ^VB ~ (VaVa) ~ (V/3Vb) + 0(e 2 )) 

= B • VB + 0(e 2 ) (45) 

If E„ = E - /U.B denotes the energy associated with motion parallel to the field 
lines, this leads to 

v c - (2cE„/eB) (b x (B • VB)) (46) 

The drifts thus derived are local and vary from point to point; however, 
making use of the constancy of /j. and E , they are easily derived if the initial 
conditions and structure of the field are given. On time scales of 0(e -1 ) gyration 
times, the accumulated drift will be of order R g , justifying eq. (38) for such time 
scales. 
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The 2nd Adiabatic Invariant J 2 

If the dependence of B on the coordinates is known, one may deduce its 
dependence on s Q and use this to integrate the motion along field lines. In par- 
ticular, if this motion turns out to be periodic, it will have an action variable 

J 2 = <^P 2 dq 2 (47) 


which may be shown from general principles^ 21 ) to be adiabatically invariant 
under slow perturbations, e.g. those arising from the drift motion. This may 
also be shown directly, as will be done in the following calculation, which re- 
sembles one by Northrop (22) but is shorter due to the canonical formulation. 

We begin by associating the particle with guiding center at 

(a, /3, s Q ) - (cpj./e, q lt q 2 ) 


with an "instantaneous adiabatic invariant" 


J 2 (E, p r q t ) 


= f { V2mE- M B /|Vs 0 | Jdq 2 
=■ f {p„/|Vs 0 |}dq 2 


(48) 


where E is the energy integral and P M is defined by the last relation. Qualita- 
tively, J 2 may be viewed as the value J 2 would assume if the particle ceased 
drifting. 
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Actually, J ' 2 will slowly change, due to its dependence on p x and q x which vary 
in the course of the drift motion. The purpose of this calculation is to show that, 
to the lowest order of expansion, this variation is periodic and averages zero. 

By equation (38) 

(Vs Q ) 2 = 1 ~ 2(a - a Q )(Va • Vs 0 ) - 2(0 - /3 Q )(Vb • Vs 0 ) 




- 1 - 2(cpj/e - a 0 )(3 a/3 s 0 ) - 2(q x - /3 0 )(3b/3 s 0 ) 

= 1 - 2u (49) 

To the' same order 

j; P„ (1 + u + (3/2)u 2 )dq 2 (50) 

Hamilton's equations give 

dj^/dt = ^(^Ji/^q,) + p 1 (3Ji/3p 1 ) (51) 

= (3H/3 Pl )(3Ji/3 qi ) - (3H/3q,)(3j'/3 Pl ) 


Explicitly 

2m(3H/3q x ) = ^(3 B/3 q x ) - 2p 2 (3 u/3 q x ) (52) 

with an analogous expression for (3H/3p x ) and similar ones for primed quanti- 
ties. Since q x and p x are associated with the drift motion, they are 0(e) and 
so are the terms of eq. (52). 

To express dj^/dt , we have to differentiate the integrand of (50), discarding 
higher order terms. In this integrand, if P M is 0(1), then u is 0(e) (see eq. 38); 
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however, taking BP n /Bqj as 0(e) (it contains a term of eq. 52), we find that 
Bu/3qj is also of the same order, due to the factor q t in eq. (49). Thus, keeping 
the two lowest orders 


= j" //BB/Bq^ 1 + u) + 2 P |T 2 (Bu/Bqj + 3i> Bu /Bq t ) J /2 P H dq 2 
= /U.3B/Bqj( 1 + u) + 2p 2 ( 1 ~ 2u) (Bu/Bqj + 3u Bu/Bqj) | /2 P„ dq 
= - m <j> j ( 1 + u) (BH/Bqj) +0(e 3 ) J /P M dq 2 (53) 


By Hamilton's equations 


q 2 = BH/Bp 2 = p 2 (l - 2u)/m 
dq 2 = P M dt/(l +u)m (54) 

If r is the period of the motion in q 2 , one may write 

= - f T (3H/3q,) dt 
Jo 

where among the arguments of H, q 2 and p 2 vary with t but q t and p x are held 
fixed. Substituting this in (51) and using Hamilton's equations 

d J 2 / d t' = J T |(BH7Sqj) (BH/Bp,) - (BH’/B Pl ) pH/Bq,) J dt (55) 

where q 2 and p 2 that appear in H are functions of t and participate in the inte- 
gration, while those in H' relate to the particle's parameters at time t 1 at which 
Bj 2 /dt is evaluated and do not participate. The expression, in general, does 
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not vanish; however, let us calculate the total change AJ 2 of J 2 in one period 

a J ' 2 = f(dj'/df) dt' 

Jo 

-= f T f’cBH'/Bqj) (3H/3 Pl ) - (BH'/3 Pl ) (BH/Bq^ dt dt' (56) 

Jo Jo 

If H were the same function of t as H' of t ' , this vanishes by symmetry, since 
an interchange of the dummy variables then reverses the sign. This, however, is 
not strictly true, for two reasons. First, the arguments q 2 (t) and p 2 (t) of H 
follow the motion of a particle the drift of which has been "turned off, whereas 
the corresponding variables in H 1 — call them q 2 (t‘) and p 2 (t') — follow the actual 
trajectory. Secondly, both H and H' depend on qj and pj of the drifting particle, 
which in turn depend on t' alone. 

However, due to the slowness of the drift, the integrand may be expanded in 
the small differences between (q t , p^ and their initial values (q 10 , p 10 ) at t = 0, 
and between (q 2 , p 2 ) and (q 2 , p 2 ). The leading term of this expansion is com- 
pletely symmetric, so that to lowest order AJ 2 vanishes. 

Let us evaluate what this means. If E is the energy, by (48) and (55) 

= O(Et) (57) 

The integrand of (56) is proportional to the energy involved in the drift motion 
and thus 0(e 2 E); however, since its lowest order does not contribute, 

Aj; = 0(e 3 Er 2 ) 
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Combining, and assuming r is of order e” 1 

AlogJ^ = 0(e 2 ) (58) 

Thus in the course of e~ l periods, one would expect J 2 to vary by no more 
than 0(e) of its value. The integral J 2 of (47) closely approximates J 2 and its 
long-term behavior is the same. 

Motion in the Earth's Magnetic Field 

An example of the preceding is afforded by radiation-belt particles moving 
in the earth's magnetic field, which may be viewed as a dipole field with external 
and internal distortions added. A typical field line starts from one hemisphere 
and ends in the other; it has a high field intensity B near its ends and a minimum 
of B somewhere near the point at which it is farthest away from earth. 

A particle trapped on such a line and conserving fi will bounce back and 
forth between "mirror points" at which 

B = = B 

m 

slowly drifting azimuthally at the same time and conserving J 2 while doing so 
(Figure 3). 

The bounce motion is evidently periodic in s, if s is defined to be measured 
from the mirror surface B(x ,y ,z) = B m , but it is not immediately evident that it 
can also be regarded as periodic in s Q , as required by the preceding formalism. 
This, however, may be ensured by a proper choice of s . By substitution in 
equations (21), it may be shown that the general transformation 
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s -> s' - s + f (a, J3) 

leads to 

s 0 ** s o = s 0 + f (a, /3) + (a - a 0 )(B f/Ba) 

+ ( r 3 “ / 5 0 )( 3 f /B ^) 

Given some function g(a,/3), one can in principle find a solution f for 

( a “ a 0 )(^ f /^a) + (fi - /. 3 0 )(B f/B/3) + f (a, /3) + g(a, /3) = 0 

which means that for a prescribed transformation of s Q 

s 0 ■* s 0 = S 0 " gC®-/ 3 ) ^ 59 ) 

there usually exists some transformation of s that leads to it. Now suppose that 
s is given as vanishing on one mirror surface. Every point on that surface can 
be labeled by a pair of values ( a,/3 ), and there will in general exist a function 

g(a,/3) giving the values assumed by s 0 on that surface at each point. If we now 

derive and apply that transformation of s which makes s 0 transform as in 

eq. (59), we will have arrived at an (a,/ 3, s 0 ) system in which s 0 vanishes on 
the mirror surface. The motion is then periodic in q 2 and all preceding 
arguments apply. 

In the geomagnetic field, the drift motion may carry the particle all the way 
around the earth — as in the case of the dipole field — or else the particle may 
drift off and eventually become "untrapped" before completing its full circuit^. 
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In the former case , the motion will have a third periodicity and an associated 
adiabatic invariant, the so-called Flux Invariant $ {6) . In the present case there 
exists no need for 4>, since m, J 2 and the energy integral E are sufficient to 
describe the motion; it comes into its own in the case of time dependent magnetic 
fields, which will generally have an associated electric field as well. It is 
possible to extend the Hamiltonian formalism to include such fields, but this is 
beyond the scope of the present work. 

Another application, which will be only briefly mentioned here, involves the 
Liouville equation corresponding to the guiding center Hamiltonian. Writing 
down this equation leads almost immediately to the gyration averaged version 
of the equation describing the behavior of a collisionless plasma, the so-called 
CGL (for Chew, Goldberger and Low) approximation to the Vlasov equation. 

Conclusion 

The preceding review of the properties and applications of Euler potentials 
is in no sense complete. The example worked out, involving the canonical 
description of guiding center motion, will provide those persistent enough to 
follow it with some practice in their application, but it is only one of a 
rapidly growing number of applications of Euler potentials. Because of this 
increasing importance, it appears that the time has come to include at least a 
brief introduction to Euler potentials in the standard course on Classical 
Electrodynamics. Hopefully, this article will provide the material and the 
background for such an introduction. 

The author is grateful to Dr. Thomas Birmingham for his critical review 
of this work. 
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SSCTIO SECUNDA DE PRINCIPII8 MOTUS FLUIDORUM 


co 


SCHOLION 2 


49. Ex solutione problematis, dum per gradus 
positam siimus progress! , casum, quo fluidi donsitas q 


acquatio ita se habet 


(£-) 


+ 


(£W£) 


o, 


ad aequationem pro- 
est constans et prima 


rcsolvcre poterimus, quod eo magis est notandum, quod huius solutionem cx 
gcncrali, quam dedimus, r’erivare iioj» licet. Quamquain autem hie ires tan turn 
variables x, y et z connderantur, tamen nihil impedit, quominus in solutions 
ilii data etiara quartum t introducamus, earn quasi conslantera spectamio. 
Sumtis ergo pro lubitu duabir ^unctionibus F , O quatuor variabilium x, y, z 
et / , ex iis ternae celeritates u, v et w ita determinabnntur, ut sit 


(dF\ (<IQ\ idF\(dG\ 
U ~ \dyj\dz) \dz)\dy) 

(dF\fdG\ (dF\[dG\ 

v -[wl{^)~[^)Kd^l 

(dF\(dG\ ( dF\ ( dG\ 

W -(^)Uy) W/W 


-f r:(y,z,t) 

-{■ A :{x, z, t) 
-f Z:{x, y, t ) , 


u'oi totum momentum iterum in eo est situm, quod cuivis membro cuiusque 
formae in reliquis respondeat aliquod, quod cum eo datum factorem habcat 
communem et signo contrario sit affectum. Totus autem hie casus, quo den- 
si tas fluidi est quantitas constans, meretur, ut seorsim diligentius evolvatur, 
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